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ABSTRACT 

We propose to use the large-scale structure of the universe as a cosmic standard ruler, based on 
the fact that the pattern of galaxy distribution should be maintained in the course of time on large 
scales. By examining the scale-dependence of the pattern in different redshift intervals it is possible 
to reconstruct the expansion history of the universe, and thus to measure the cosmological parameters 
governing the expansion of the universe. The features in the galaxy distribution that can be used as 
standard rulers include the topology of large-scale structure and the overall shapes of galaxy power 
spectrum and correlation function. The genus, being an intrinsic topology measure, is resistant against 
the non-linear gravitational evolution, galaxy biasing, and redshift-space distortion effects, and thus 
is ideal for quantifying the primordial topology of the large-scale structure. The expansion history of 
the universe can be constrained by comparing among the genus measured at different redshifts. In 
the case of initially Gaussian fluctuations the genus accurately recovers the slope of the primordial 
power spectrum near the smoothing scale, and the expansion history can be constrained by comparing 
between the predicted and measured genus. 

Subject headings: large-scale structure of the universe - cosmology: theory 



L INTRODUCTION 

There are three kinds of phenomena of the universe 
that are currently used to constrain cosmological mod- 
els. The first is the primordial fluctuations or the initial 
conditions. Currently available tracers of the primor- 
dial fluctuations are the cosmic microwave background 
(hereafter CMB) anisotropies and the large-scale struc- 
ture (LSS) of the universe. From these one can study the 
geometry of space, matter content, matter power spec- 
trum (PS), non-Gaussianity of the initial conditions, and 
so on. Information from these tracers has limitations 
because it contains knowledge only in one thin shell lo- 
cated at a specific epoch in the case of CMB, or because 
the amount of the corresponding data is not yet large 
enough to constrain cosmological models strongly in the 
LSS case. The eventual limitation lies in the finite vol- 
ume of the observable universe. 

The second measurable phenomenon of the universe is 
the expansion of the universe. It can be measured by 
observing the standard candles (e.g. supernova type la; 
Colgate 1979; Riess et al. 1998; Permultter et al. 1999), 
the standard rulers (e.g. baryon acoustic oscillations, 
hereafter BAOs; Peebles & Yu 1970; Meiksin, White, & 
Peacock 1999), or standard populations, if any. Red- 
shifts of these objects give us the relation between the 
comoving distance r and redshift z through the lumi- 
nosity distance Dl{z) and/or angular diameter distance 
Da{z)^ which constrain the expansion history of the uni- 
verse or the Hubble parameter H{z) through the relation 
r{z) = /q dz'/H{z'). The Hubble parameter depends on 
many cosmological parameters such as the total density 
parameter f2totj matter density parameter Qm, and the 
equation of state of the dark energy w = P/p. However, 
there are various kinds of systematic effects that limit the 
power of this method. For example, the dependences of 
the 'standard' properties on tracer subclasses and on red- 
Electronic address: |cbp@kias.re.kr, yrk@kias.re.kr| 



shift are the most serious error sources in measuring r(z) 
in the case of the standard candles and populations. The 
standard rulers also suffer from all kinds of systematics 
such as non-linear gravitational evolution, redshift-space 
distortion, past light-cone effects, and biasing of tracers. 

The third phenomenon is the growth of cosmic struc- 
tures, which depends on both expansion history and ini- 
tial matter fluctuations. This can be examined by ob- 
serving the integrated Sachs- Wolfe effect causing a corre- 
lation between CMB anisotropy and LSS (Sachs & Wolfe 
1967, Corasaniti et al. 2003), abundance of galaxy clus- 
ters (Allen et al. 2004, Rapetti et al. 2005), and the 
weak gravitational lensing by LSS (Cooray & Huterer 
1999). Properties of some non-linear objects can be also 
used. Various present and redshift-dependent properties 
of intergalactic medium (near the reionization epoch, in 
particular), massive dark halos (luminous galaxies and 
clusters of galaxies), etc., are the combined results of the 
initial matter fluctuations, expansion history, and non- 
linear physics. 

In this paper we propose to use the pattern of the 
large-scale galaxy distribution to study both the first 
and second phenomena of the universe. We will intro- 
duce this tool as a geometrical method similar to the 
Alcock-Paczynski test (Alcock & Paczynski 1979) or the 
B AO-scale method (Blake & Glazebrook 2003). In the 
forthcoming papers we will also show that this method 
is complementary to other methods such as the BAO- 
scale method, and has a power comparable to the BAO 
method in constraining the dark energy equation of state 
(Kim et al. 2009). 

2. LARGE-SCALE STRUCTURE AS A STANDARD 

RULER 

The large-scale distribution of galaxies has long been 
used to constrain cosmological models through two-point 
correlation function (hereafter CF; Davis & Peebles 1983; 
Maddox et al. 1990) and PS analyses (Park, Gott, & 
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da Costa 1992; Vogclcy ct al. 1992; Park ct al. 1994, 
Tegmark ct al. 2006) because the shapes of the PS 
and CF depend on the cosmological parameters such as 
the matter and baryon density parameters {flmh'^, fifo/i^), 
and the primordial spectral index (rig). 

Sensitivity to the expansion of the universe appears 
when the LSS observation spans a range of redshift. 
Comparing of the shape of the PS at two different epochs, 
knowing they should be the same, one can find how 
the universe has expanded between the epochs. This 
would have been impossible if the universe had a scale- 
free PS since there is no characteristic scale to compare 
the spectra at two epochs. But the curvature (the scale- 
dependence of the slope) in the observed CF (Maddox 
et al. 1990) and the PS (Vogeley et al. 1992; Park et 
al. 1994) has been well confirmed. Theoretically, the 
linear density PS of the Cold Dark Matter models is ex- 
pected to have a peak near the scale corresponding to 
the epoch of matter-radiation equality, approaches fc"^ 
at the largest scales, and at the smallest scales. 

The PS or CF is only one of many properties of LSS 
that can be used as standard rulers, and another is the 
topology. Bond et al. (1996) pointed out that the 
filament-dominated cosmic web is present in embryonic 
form in the overdensity pattern of the initial fluctua- 
tions with non-linear dynamics just sharpening the im- 
age. However, it is not just overdensity pattern like the 
cosmic web structure hut the whole large-scale structure 
including the cosmic voids that mcinorizes their initial 
birth places. After all, it is the whole cosmic sponge 
rather than just its high density side that remains un- 
changed through the evolution of the universe. In the 
case of the flat ACDM model with the WMAP 5-year 
cosmological parameters of ^lm = 0.26, i^A = 0.74, ili, = 
0.044, /i = 0.72,0-8 = 0.796, and = 0.96 (Dunkley et 
al. 2009), our N-body simulations show that the RMS 
matter displacements till redshifts z = 0.5 and are 
7.7 and 9.7 ft.~^Mpc, respectively. At the scales much 
larger than the RMS displacement the topology of LSS 
should represent that of the initial density fluctuations 
accurately. 

In this paper we will adopt the genus statistic (Gott, 
Dickinson & Melott 1986) as a measure of the LSS topol- 
ogy. For Gaussian random phase initial conditions the 
genus curve is given by 

g{u)^A{l-i.^)e-'''/^ (1) 

where ly is the density threshold level normalized by 
the RMS density fluctuation, and the amplitude A = 
((A;^)/3)^/^/27r^ and (A:^) is the average value of /c^ in the 
smooth PS (Hamilton, Gott & Weinberg 1986; Doroshke- 
vich 1970). The amplitude A measures the slope of the 
PS near the smoothing scale, and is independent of the 
amplitude of the PS. According to most previous studies 
of non-Gaussianity of the primordial density fluctuations, 
the large-scale distribution of galaxies is consistent with 
the cosmological models with initially-Gaussian matter 
density fields (e.g. Gott et al. 2009). 

The topology of LSS at a given scale is conserved in 
the course of time in the comoving space regardless of 
whether the initial topology is Gaussian or non-Gaussian 
as long as the scale corresponds to the linear regime. 
When the primordial field is not Gaussian, A is not sim- 
ply determined by the shape of the PS. However, A is still 
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Fig. 1. — Amplitude of the genus curve per unit smoothing 
volume for the flat ACDM models with Qm + f^A = 1- The label 
on each line is Qrn and the remaining cosmological parameters are 
set to those from the WMAP 5-year data. The thick line is for 
f!m = 0.26, and the boundaries of the shaded region correspond to 
f^m = 0.29 and 0.23, which are chosen from the mean and the 1<t 
uncertainty limits of the WMAP 5-year data. The data points are 
from the SDSS DR4plus sample (Gott et al. 2009). 

a conscirvcid quantity and can be used for reconstructing 
the expansion history of the universe. In this sense the 
topology of LSS is a cosmic standard ruler independent 
of the PS. 

The genus is measured from the iso-density contour 

stirfaces of the smoothed galaxy distribution. Being a 
measure of intrinsic topology, the genus is insensitive to 
the galax;y biasing and rcdsliift-space distortions. This 
is because the intrinsic topology does not change as 
the shape of a structure continuously deforms without 
breaking up or connecting with itself or other structures. 
Furthermore, according to the second-order perturbation 
theory, there is no change in A due to the weak non-linear 
gravitational evolution (Matsubara 1994). Namely, the 
genus amplitude is a powerful measure of the slope of 
the primordial PS in the case of a Gaussian field. It is 
because of this property of the genus topology that we 
prefer to adopt it as a standard ruler rather than directly 
using the PS or CF. 

Figure 1 shows the amplitude of the genus curve per 
unit smoothing volume, gRQ, when the matter field is 
smoothed over Rq by a Gaussian filter, for a series of the 
flat ACDM models with O^ + Oa = 1 and /i = 0.72. The 
label on each curve is the value of Qm- Each curve has a 
characteristic shape reflecting the shape of the ACDM PS 
with different f^m- A comparison between these curves 
and gRQ measured from low redshift LSS data constrains 
the cosmological parameters related with the shape of the 
PS. The uncertainty in from the WMAP 5-year data 
corresponds to 6.6% variation in the genus amplitude at 
Rg = 20/i~^Mpc, for example (see the shaded region in 
Fig. 1). The two data points are from Gott et al. (2009) 
who mcastirecl A with 4% uncertainty at 21 h~^Mpc us- 
ing the SDSS DR4plus sample. 

3. EXPANSION HISTORY OF THE UNIVERSE 
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Fig. 2. — A schematic diagram illustrating changes in the unit 
volume (dashed boxes) and smoothing scale (error bars) when dif- 
ferent r{z) relations are adopted for large-scale structure data ob- 
served at low and high redshifts, supposing the true cosmology is 
X and the assumed cosmologies are Y and Z with wrong w's. 



As the universe expands, the cosmic sponge remains 
unchanged and the shape of the PS at large scales is con- 
served in the comoving space. If an observational sam- 
ple of LSS covers a range of redshift, one can measure 
the PS or CF in different redshift intervals and compare 
their shapes to check whether or not the adopted r{z) 
relation is correct. If their shapes in different redshift 
intervals do not agree with one another, we change the 
cosmological parameters and thus r{z) relation until an 
agreement is achieved. In this section we will consider 
only the Gaussian fluctuation case, where the genus am- 
plitude gives information equivalent to the shape of the 
PS, for a demonstration of concept. 

Figure 2 illustrates what happens to the LSS analy- 
sis when wrong r{z) relations (wrong cosmologies) are 
adopted. Suppose one is trying to constrain w while 
other parameters are fixed. The box on the left shows 
LSS at low redshift in a square region of the universe 
with w = —I. The error bar indicates the smoothing 
scale. If one transforms the redshifts of distant galaxies 
to comoving distances using the r{z) relation adopting 
the correct value of w, a region having the same comov- 
ing size located at high redshift will enclose the same 
comoving volume and the smoothing length will corre- 
spond to the scale equal to that at z = (the middle box 
on the right). But if we choose w = —0.5, we mistakenly 
think the space is expanding slower than the reality. As 
a result of the wrong r-z transformation, a unit volume 
(the dashed box in Fig. 2) will enclose more LSS than 
the box does at 2 = 0, and the smoothing length cor- 
responds to a scale larger than that at low z. Because 
the box of a unit comoving volume contains more LSS 
but the smoothing is made over a larger scale, their ef- 
fects on the genus partially cancel with each other but 
there remains some net cffcict unless the density fluctua- 
tion field is scale-free. When one adopts a universe that 
expands faster than the real one (w = —1.5 for example), 
the comoving volume of the box at high redshift actually 
amounts to a smaller volume compared to that in the 
true cosmology, and the smoothing scale corresponds to 



a scale smaller than what it is intended to. 

The amplitude of the genus curve when a wrong cos- 
mology 'Y' is adopted while the true cosmology is 'X', 
can be estimated in the following way. The volume factor 
at redshift ^; in a cosmology is given hy V = D\/H{z) 
(Peebles 1993). When a wrong cosmology is adopted, 
the fractional change in volume is Vx{z)/Vy{z) when 
the samples are constrained to have the same comov- 
ing volume under the given cosmologies. (Due to the 
Alcock-Paczynski effect, the amormt of radial and tan- 
gential length variations is slightly different from each 
other. In the present treatment we average structures 
over angles and consider only the volume effect.) On the 
other hand, the smoothing length changes by a factor 
Xxviz) = [Vx(^)/Vy(z)]i/3. Therefore, the amplitude 
of the genus curve measured when the wrong cosmology 
Y is adopted to convert redshifts to comoving distance, 
becomes 



gyiz; Rg) = 9x{z; R'G)Vx{z)/Vy{z), 



(2) 



where R'q = \xy{z)Rg- This formula is equivalent to 

gyiz; Rg)Rg = 9x{z: R'q)R'%. Any change in cosmology 
that affects the expansion history of the universe will re- 
sult in change in the redshift dependence of g{RG)- One 
can use this formula to estimate the genus in a partic- 
ular cosmology when its value for a fiducial cosmology 
is known without making full numerical mock sample 
analysis. Cosmological parameters can be constrained 
through an iterative process to minimize the difference 
between observations and the theoretical prediction. In 
practice, observational data occupy a finite redshift in- 
terval, and the scaling factors in Equation (2) should 
be replaced by an integral over redshift. It should also 
be pointed out that the number density of the objects 
tracing the LSS will show a radial gradient due to the 
redshift dependence of the volume factor when a wrong 
cosmology is adopted. If the radial distribution of the 
tracers is forced to be uniform, the selection criterion of 
the objects will become non-uniform. 

When one measures the two-point CF under the as- 
sumption of a wrong cosmology y, one will incorrectly 
scale the separation between galaxies, and the posi- 
tion of features in the CF moves by the scaling fac- 
tor, namely ^y{s]z) = ^x{s';z) where s' = \xy{z)s 
(see also Percival et al. 2007). This causes the slope 
of the CF to change. Likewise, the PS is scaled as 
Pr(fc;z) = Px{k';z) where k' = X^\.{z)k. 

Figure 3 shows how the genus is used to measure the 
cosmological parameters that govern the expansion his- 
tory of the universe. Here it is again assumed that the 
true cosmology has 0,^ = 0.26 and Oe = 0.74 with 
w = — 1 (The subscript E stands for the dark energy). 
We adopt Qm + — 1, and all remaining parameters 
are fixed to the WMAP 5-year parameters. The left 
panel shows the predicted amplitude of the genus curve 
at Rg = 15 h~^Mpc when different sets of (Clin,w) are 
adopted. Each line is independent of z because the shape 
of the linear PS is conserved in the comoving space. It 
is also independent of w since w does not change the 
shape of the linear PS. The lines in the right panel are 
the genus per unit smoothing volume that will be mea- 
sured when five different cosmologies are adopted. It 
demonstrates that the redshift-dependence of the genus 
amplitude is quite different for cosmologies with differ- 
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Fig. 3, — The genus per unit smoothing volume at Rq = 
15/i~^Mpc for five flat cosmological models with different sets of 
(n^jUi) (left panel). The right panel shows the genus that will be 
actually measured when five different sets are assumed in the r-z 
transformation even though the correct one is (0,26,-1,0). 

ent rim and w. The correct choice of (0.26, —1.0) will 
result in an agreement between the theoretical predic- 
tion, the horizontal line with gR^ = 0.00470, and the 
obscrvationally measured one. But when (0.26, —0.5) 
is mistakenly adopted, the wrong r{z) makes the genus 
amplitude overestimated. This means the volume factor 
dominates the change in Equation 2. At z = 0.5, the 
measured value will be 5.3% larger than the predicted 
value. When the fractional error in the observed g is Ag, 
the fractional error in w constrained by this single data 
point is roughly A„ = 9.5Ag. Therefore, if the genus 
is measured with accuracy better than 1% at z ^ 0.5, 
one can constrain w with error less than about 10% by 
comparing the theoretical prediction with what is actu- 
ally measured. When a model with more negative w is 
adopted, the measured genus amplitude falls below the 
predicted value [see the (0.26, —1.5) case]. 



4. SUMMARY 

The key points of this paper can be summarized as 
follows: 



1. We point out that the pattern of the LSS is con- 
served in time and can be used as a cosmic standard 
ruler. 

2. We propose to use the sponge topology of LSS as a 
more robust standard ruler than P{k) or ^(r) be- 
cause the intrinsic topology statistics are less sus- 
ceptible to various non-linear systematic effects. 

To enhance the power of the standard rulers it is nec- 
essary to have accurate knowledge on the systematic ef- 
fects such as non-linear gravitational evolution, scale- 
dependent galaxy biasing, and redshift-space distortion 
(Park, Kim & Gott 2005; James et al. 2009). This will 
be the subject of our subsequent paper (Kim et al. 2009). 

Recently, Gott ct al. (2009) measured the genus ampli- 
tude with 4% error at the smoothing scale of 21 /i~^Mpc 
from the Luminous Red Galaxy sample of the SDSS 
DR4plus. The current redshift survey sample already 
reached the size that enables us to do the topology study 
at a few percent uncertainty level. The final SDSS DR7 
data and the future LSS surveys are expected to signif- 
icantly increase the accuracy. Various constrains on w 
from existing and future surveys will be presented by 
Kim et al. (2009) 

In addition to the 3D genus, one can also use the scale- 
dependence of the ID level-crossing statistic or the 2D 
genus of LSS as the standard rulers. For example, cos- 
mological parameters can be estimated by requiring the 
level crossings per unit comoving length in the radial 
distribution of Ly-a forest clouds to be constant in time. 
The 2D genus of the galaxy distribution in the photo- 
metric redshift slices can be also used to constrain the 
expansion history of the universe. We will explore the 
usefulness of these statistics in the forthcoming papers. 
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